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Abstract. We consider discrete and integer-valued cellular automata (CA). A particular 
class of which comprises “Hamiltonian CA” with equations of motion that bear similarities 
to Hamilton’s equations, while they present discrete updating rules. The dynamics is linear , 
quite similar to unitary evolution described by the Schrodinger equation. This has been 
essential in our construction of an invertible map between such CA and continuous quantum 
mechanical models, which incorporate a fundamental discreteness scale. Based on Shannon’s 
sampling theory, it leads, for example, to a one-to-one relation between quantum mechanical 
and CA conservation laws. The important issue of linearity of the theory is examined here by 
incorporating higher-order nonlinearities into the underlying action. These produce inconsistent 
nonlocal (in time) effects when trying to describe continuously such nonlinear CA. Therefore, in 
the present framework, only linear CA and local quantum mechanical dynamics are compatible. 


1. Introduction 

A novel analysis of quantum mechanics (QM) — which aims at redesigning the foundations of 
quantum theory — has recently been proposed by G. ’t Hooft [lj. The hope for a comprehensive 
theory expressed in this far-reaching considerations is founded on the observation that quantum 
mechanical features arise in a large variety of deterministic “mechanical” models. While 
practically all of these models have been singular cases, i.e., which cannot easily be generalized 
to cover a realistic range of phenomena incorporating interactions, CA promise to provide the 
necessary versatility mm®. For an incomplete list of various earlier attempts in this field, see, 
for example, Refs. ISU6UTU8U91 nni hd na [13 and further references therein. 

The linearity of quantum mechanics (QM) is a fundamental feature, which is particularly 
visible in the Schrodinger equation (leaving aside here models which attempt to describe 
measurement processes dynamically). This is independent of the particular object under 
study, provided it is sufficiently isolated from anything else, and is naturally reflected in 
the superposition principle. Thus, linearity entails the “quantum essentials” interference and 
entanglement. 

The linearity of QM has been questioned and nonlinear modifications have been proposed 
earlier — not only as suitable approximations for complicated many-body dynamics, but 
especially in order to test experimentally the robustness of QM against such nonlinear 
deformations. This has been thoroughly discussed by T.F. Jordan who presented a proof ‘from 
within’ quantum theory that the theory has to be linear, given the separability assumption 


that the system we are considering can be described as part of a larger system without 
interaction with the rest of the larger system.” |14j 

In our recent work, we have considered a discrete dynamical theory , which deviates notably 
from quantum theory, at first sight. However, we have shown with the help of sampling theory 
that the deterministic mechanics of the class of Hamiltonian CA can be related to QM in the 
presence of a fundamental time scale. This relation appears to demonstrate that consistency of 
the action principle of the underlying discrete dynamics implies, in particular, the linearity of 
both theories. 

We will review some pertinent results, followed by a more detailed study of the undesirable 
consequences that we are confronted with, if we generalize the action principle by incorporating 
genuine nonlinearities and try to enlarge the proposed class of CA in this way. 

The CA approach may offer additional insight into interference and entanglement, in the 
limit where the discreteness scale can be treated as sufficiently small. Furthermore, future 
developments are conceivable which can address the dynamics of QM measurement processes, 
about which standard quantum theory remains silent. 


2. An action principle for cellular automata 

In the following, we recall briefly the dynamics of a class of discrete CA automata introduced 
earlier, which bears remarkable similarity with Hamiltonian dynamics in the continuum on one 
side and with QM on the other 

We describe the state of a classical CA with countably many degrees of freedom by discrete 
integer-valued “coordinates” x“,r n and “conjugated momenta” p“,7r n , where a £ No denote 
different degrees of freedom and n £ Z different states. The x n and p n might be higher¬ 
dimensional vectors, while r n and V n are assumed one-dimensional. — The “coordinate” r n has 
been separated from the x“’s (correspondingly ir n from the s), since this degree of freedom 
allows to represent a dynamical time variable, discussed in Enanansi, with further references 
given there. 

For any one of the dynamical variables, say f n , the finite difference operator A is defined by: 

A/n := fn - fn—l ■ (1) 

Furthermore, we introduce the quantities (assuming the summation convention for Greek indices, 


tin •— C n 7T n , (2) 

H n := ^Sapip^ + xZxfo+AappZxP + Rn , (3) 

A n := AT n (i? n + H n ~ i) + a n , (4) 

where the constants, c n , the symmetric, S = {S^}, and the antisymmetric, A = {A a p}, 


matrices are all integer-valued. The definition @ determines the behaviour of the variable r n 
and only the most simple choice (involving a single constant) will be relevant for our purposes, 
cf. below. R n stands for higher than second powers in or p". — Discussion of such genuine 
nonlinearities is the aim of the present note and will be resumed in due course. 

Based on these definitions, we introduce the integer-valued CA action by: 

<S := + Pn-l) Ax n + fa* + 7T n _i)A T n - An] . (5) 

n 

Furthermore, we consider integer-valued variations 5f n to be applied to a polynomial g: 

^fridifn) ■■= [g(fn + & fn) ~ htfn ~ 5f n )\/25f n , 


( 6 ) 


and 5f n g = 0, if 5f n = 0. In terms of these notions, we postulate the following variational 
principle for the CA dynamics. 

The CA Action Principle: The discrete evolution of a CA is determined by the stationarity of 
its action under arbitrary integer-valued variations of all dynamical variables, SS = 0. • 

Some characteristic features of this CA Action Principle deserve to be mentioned: 

i) Variations of terms that are constant, linear, or quadratic in integer-valued variables yield 
analogous results as standard infinitesimal variations of corresponding terms in the continuum. 

ii) While infinitesimal variations do not conform with integer valuedness, there is a priori no 
restriction of integer variations, hence arbitrary integer-valued variations must be admitted. 

iii) However, for arbitrary variations 5f n , the remainder of higher powers , R n in H n , which 
enters the action, has to vanish for consistency. Otherwise the number of equations of motion 
generated by variation of the action, according to Eq.(j6|). would exceed the number of variables. 
(A suitably chosen Rq or a sufficient small number of such remainder terms can serve to encode 
the initial conditions for CA dynamics.) 

As we observed in earlier work, these features seem to be essential when constructing a 
map between the considered CA and equivalent quantum mechanical models based on sampling 
theory, cf. Sec. 3. We will try to illuminate this in Sec. 4., by studying a generalization of 
the variational principle incorporating nonlinearities and by pointing out consequences of such 
nonlinearities which obstruct this way to arrive at QM models from CA. 

2.1. The CA equations of motion and conservation laws 

We now apply the CA Action Principle to the action S (keeping R n = 0 for the moment), which 
yields the following discrete equations of motion: 

®n = Tn(S a/3 pP + A a pX^) , = -T n (S al3 xP-A al3 pP) , (7) 

l~n — C n , 7 T n — H n , ( 8 ) 

introducing the suggestive notation O n := O n+ \ — O n -\ . We emphasize that all terms herein 
are integer-valued. The fact that we arrive at finite difference equations reflects the discreteness 
of the automaton time n and their appearance has motivated the name Hamiltonian CA. 

The equations of motion are time reversal invariant , since the state n + 1 can be calculated 
from knowledge of the earlier states n and n — 1 and the state n — 1 from the later ones n +1 and 
n. — Note that the f n present parameters for the evolving x,p- variables, as a consequence of 
Eqs. (jSJ). More generally, t can play the role of a dynamically coupled lapse function in Eqs. (0. 

Introducing the self-adjoint matrix H := S + iA, the Eqs. 0 can be combined into: 

Xn + iPn = -if n H a p{x^ n + ip?) , (9) 

and its adjoint. Thus, we find here a discrete analogue of Schrodinger’s equation , with 
:= + ip^ as amplitude of the “a-component” of “state vector” | if) at “time” n and 

with H as Hamilton operator. (We will use QM terminology freely here and in the following.) 

Furthermore, there are conservation laws respected by the discrete equations of motion, or 
by Eq. Q, which are in one-to-one correspondence with those of the corresponding Schrodinger 
equation in the continuum In particular, the Eqs. 0 imply the following theorem. 

Theorem A : For any matrix G that commutes with H , [G,H] = 0, there is a discrete 
conservation law: 

r n a G a ^ n +r n a G a ^ n = o. 


(10) 




For self-adjoint G, with complex integer elements, this relation concerns real integers. • 

Corollary A: For G := 1, the Eq. (11011 implies a conserved constraint on the state variables: 

O“ + On =0 • (H) 

For G \= H , an energy conservation law follows. • 

Note that Ec|s. (j 10 [) and m cannot be “integrated” as usual, since the Leibniz rule is 
modified here. Recalling O n := O n+ \ — O n -i, we find, for example, O n+ iO' n+1 — O n -iO' n _ l = 
lj(O n [0' n+1 + O r n _ i] + [O n+ 1 + 0 n -i\0' n ), instead of the product rule of differentiation. 

Furthermore, the continuum limit of the equations of motion and their conservation laws 
does not simply follow from letting the discreteness scale l —>■ 0. The integer valuedness of 
all quantities conflicts with continuous time and related derivatives. Hence, we need a more 
elaborate mapping, in order to relate CA to continuum models. 

In the following Sec. 3., it will be shown that such an invertible map between the descriptions 
of discrete time Hamiltonian CA and of quantum mechanical objects evolving in continuous 
time can indeed be constructed, taking into account a fundamental discreteness scale l. 

3. The CA QM map based on sampling theory 

Despite the notable similarities of the Hamiltonian CA with QM systems, we may wonder 
whether the discreteness of a CA can be reconciled with a continuum description at all and, in 
particular, with QM? 

We have suggested earlier that especially wave functions, like other fields, could be 
simultaneously discrete and continuous, represented by sufficiently smooth functions containing 
a finite density of degrees of freedom mm- Related ideas have been discussed by T.D. Lee and 
collaborators and recently by A. Kempf in attempts to introduce a covariant ultraviolet cut-off 
into quantum field theories and, last not least, for gravity, see, for example, Refs. Dana with 
further references there (and in [2j[3]). However, integer-valued CA have first been connected 
with structure of QM from this perspective in our work. 

Information can have simultaneously continuous and discrete character as pointed out by C.E. 
Shannon in his pioneering work [18]. This is routinely applied in signal processing, converting 
analog to digital encoding and vice versa. Sampling theory demonstrates that a bandlimited 
signal can be perfectly reconstructed, provided discrete samples of it are taken at the rate of at 
least twice the band limit (Nyquist rate) — for an extensive review, see Ref. [19] . 

We shall make use of the basic version m of the Sampling Theorem: 

Consider square integrable bandlimited functions f, i.e., which can be represented as f(t) = 
(27r) _1 f dw e~ lult f(uj), with bandwidth uj max - Given the set of amplitudes {f(t n )} for the 

set {t n } of equidistantly spaced times (spacing TT/ui max ), the function / is obtained for all t by: 

/(*) = £ f(tn) ‘ y 1 . (12) 

n Wmax (f t n ) 

Since the CA state is labelled by the integer n, the automaton time, the corresponding 
discrete physical time is obtained by multiplying with the fundamental scale l, t n = nl, and the 
bandwidth by u max = 7r/Z. 

When attempting to map invertibly Eqs. (]7|) on continuum equations, according to Eq. (fT2lh 
the nonlinearity on the right-hand sides is problematic, since the product of two functions, with 
bandwidth u) max each, is not a function with the same bandwidth. Therefore, we assume here 
that T n is a constant and postpone the discussion of generic nonlinearities to Sec. 4. 





Recalling the discrete time equation Q, we insert iff := xf + ipf, as before, and apply the 
Sampling Theorem to obtain the equivalent continuous time equation : 

(£>i - D-i)ip a (t) = 2sinh (ld t )if a (t) = -H a ^(t) , (13) 

i 

employing the translation operator defined by Dxfit) '■= f(t + T ) and implementing the natural 
choice f n = 1, for all n. 

It appears that we recover the Schrodinger equation. However, it is modified in important 
ways, reflecting the presence of the scale l. 

First of all, by construction, the continuous time wave function if a is bandlimited (by uimax)- 
Therefore, knowing the wave function at the discrete times of a set {to + nl\n £ Z}, with to 
arbitrary, it can be reconstructed for all times by a slight generalization of Eq. (1121) . Furthermore, 
since Ea. (fT3|) is of the form f(t +l ) = f(t — l) — iHf(t), it is sufficient to know / at two times, 
say to and to — l, in order to obtain it for all times of the set {to + nl\n £ Z }. Thus, we learn 
that two initial conditions (separated by a time step l) have to be specified to define the solution 
of Eq. (1131) . This agrees precisely with the requirements of the discrete description of the CA, 
cf. Sec. 2.1. — In order that sampling reproduces the integer-valuedness of the CA, both initial 
values have to be integer-valued. 

If instead the modified Schrodinger equation is written in terms of the infinite series of odd 
powers of time derivatives, it might give the false impression that an infinity of initial conditions 
are required. In any case, the higher-order derivatives are negligible for low-energy wave 
functions, which vary little with respect to the cut-off scale, i.e. \d k if/dt k \ <C l~ k = {cvmax/^) k - 

Secondly, the bandlimit |w| < uj max leads to an ultraviolet cut-off of the energy E of stationary 
states of the generic form ifE(t) := exp(—iEt)if. Diagonalizing the self-adjoint Hamiltonian, 
yields an eigenvalue equation and a modified dispersion relation : 

E a = r 1 arcsin(e Q /2) = (2l)~ 1 e a [l + eJ/24 + 0(e*)] , (14) 

where {e Q } is the set of eigenvalues of the Hamiltonian; e.g., a could label the momentum modes 
of a given spatial lattice. We find that the spectrum {E a } is cut off by \E a \ < tt/ 21 = u max /2, 
i.e. half the bandlimit. We have discussed further aspects of this result elsewhere mm- 

Finally, the discrete CA conservation laws, Theorem A and Corollary A in Sec. 2.1., Eqs. m 
and (ED, respectively, naturally have a counterpart in the continuum description obtained with 
the help of sampling theory. The absence of an ordinary time derivative in the discrete CA model, 
where only finite differences can play a role, leads to similar combinations of translation operators 
as on the left-hand side of the modified Schrodinger equation, Eq. (1131) . in the continuous time 
conservation laws. See Refs. mm, where also related symmetries have been addressed. 

4. The options for nonlinear Hamiltonian cellular automata 

The general properties and certain quantum features, in particular, of the Hamiltonian CA that 
we recalled in the previous sections derive from the CA Action Principle introduced in Sec. 2. 
A most important aspect that we observed has been the linearity of the equations of motion, 
Eqs. © or Eq. @. 

We have repeatedly pointed out that additional higher-order terms in the action, which would 
lead to nonlinear (in xf, pf or iff) terms in the equations of motion, would simultaneously enlarge 
the set of equations. This differs markedly from what one is used to in applications of variational 
principles based on the continuum, where ordinary differential calculus is available. Presently, 
due to the absence of infinitesimals and ensuing necessity to admit arbitrary variations, the CA 
dynamics tends to become overdetermined. — As we shall see in the following, however, a closer 
look at this problem reveals some interesting information on the availability and consequences 
of nonlinear extensions. 


4-1- Generalizing the variational derivative 

The additonal equations of motion, which can render the equations of motion inconsistent since 
overdetermined, are related to two aspects of the dynamics. One comprises the additional higher 
than quadratic powers of dynamical variables in the action, which kind of terms we summarized 
by R n in definition ([3]). While the other consists in the arbitrary integer-valued variations 5f n 
of all dynamical variables present in the action, which are allowed by the variational principle. 
They are applied according to the definition of Eq. ([U]) and this can produce additional terms 
which involve powers of 5f n . The coefficients of such terms all have to vanish independently, 
leading to additonal equations of motion. The resulting enlarged set of equations, as compared 
to a given number of variables, will be overdetermined in generic cases. 

This problem can be avoided by suitably generalizing the definition of the variations. 
Replacing Eq. 0, we define here: 

Sf9 {N \f) ■= ^lk[g {N) {f + m k 5f) - g {N \f -m k Sf)\/2Sf , (15) 

fc>i 


where / stands for a dynamical variable entering the IV-th degree polynomial g( N ^ and 'y k 
and rri k (m k 7 ^ m k >, for k 7 ^ k') are constant real and positive integer-valued coefficients, 
respectively, to be determined. — Namely, we aim to arrange these coefficients in such a way 
that S f g {N) (f) = g( N 1 )(/); thus, the variation results in a polynomial g( N ^ of degree N — 1 
and all other possible terms proportional to powers of Sf cancel by construction. 

This eliminates the possibility of having an overdetermined set of (generally nonlinear) 
equations of motion. — In order to proceed, we write the polynomial explicitly, 

9 {N \f) ■= 9o + S 1/ 1 + ••• + 9n f N , (16) 


and expand the difference appearing in Eq. (1151) : 


\g {N \f + m k 5f) 


~ 9 {N \f ~ m k 5f)\/2 = gi-m k Sf + g2-^rn k fSf 
+ 93 ■ (Sm k f 2 Sf + (m fc < 5 /) 3 ) 

+ 9a ■ (4m k f 3 Sf + 4f(m k 6ff) 

+ g 5 ■ (5 m k f A bf + 10 f 2 {m k Sf) 3 + ( m k 5ff^j 
+ .... 


(17) 


Note that the terms oc Sf correspond to the ones obtained by ordinary differentiation of the 
polynomial. However, since any one dynamical variable / of our CA is integer-valued, such 
derivatives are to be interpreted only as a symbolical notation in the following. Thus, we have 
5f9 {N \f) = 12k • (d/d f)g^ N \f) + ... , where the additional terms involving powers of Sf 

are not spelled out. The point of our considerations is that the latter can be made to vanish 
always by suitably adjusting the coefficients 7 *, and m k , if we restrict the maximal order of 
polynomials to be dealt with. 

For illustration, we consider all polynomials of order < 4, i.e. gG). Here, the terms oc Sf 3 
are eliminated always, cf. Eq. (|17|) . if the following condition is fullfilled: 

^Zlk{m k f = 0 , 

k>\ 


( 18 ) 


A solution is provided by: mi = 1, m 2 = m > 2, 71 = 1/(1 — m~ 2 ), 72 = —nW 3 /( 1 — m -2 ), and 
all other coefficients vanishing. Thus, we obtain: (/) = (d/d/)^^ 4 ^ (/), c/. Ea. (fT5j) . This 

solution is sufficient for our purposes but not unique. 

By this elementary reasoning, we have obtained a satisfactory variational derivative, which 
avoids the problem of arriving at an overdetermined set of equations of motion. Our approach 
can be generalized to polynomials of arbitrary finite order; limitation to finite order being 
warranted by integer-valuedness of the variables. 

Consequently, a suitably generalized variational derivative can be employed in the CA Action 
Principle , such that consistent finite difference equations of motion incorporating nonlinear 
potentials result, which maintain the structural similarity with classical Hamilton’s equations 
that we have seen in Sec. 2 . 1 . 

4.2. Problems that arise with nonlinearity 

The generalized variational derivative of Eq. (|15l) can be employed to define a Poisson bracket , 
similarly as we discussed elsewhere ® 12 Q] — Since this variational derivative acts practically like 
an ordinary derivative, there is apparently no need to restrict the related algebra of observables 
in any way. 

Yet it does remain consistent to consider only (linear or) quadratic forms in the dynamical 
variables, recovering the previous results and analogous symmetry properties as in QM [201 EQ. 

However, once higher order polynomials (in or if®) are admitted in the action and 

equations of motion, or as observables, it is not consistent to limit the set of relevant polynomials 
at any finite order. For example, the Poisson bracket of two polynomials of order N and N', 
respectively, may result in a polynomial of order N + N' — 2 > N,N'. This is problematic, since 
arbitrarily high powers of integer-valued variables, and linear combinations thereof, generated 
in this way, will eventually lead to divergent quantitiesl 

Thus, we observe here a qualitatively profound ‘bifurcation’ in the properties of Hamiltonian 
CA tied to the presence or absence of nonlinearities in their equations of motion. The preceding 
remarks seem to imply that an algebra of observables cannot even be defined for the nonlinear 
case in an analogous way as in classical mechanics. To emphasize this point, we note that the 
kind of discrete or continuous conservation laws (and traces of QM unitary symmetry ), discussed 
in Secs. 2.1. or 3., see also Ref. M, will generally be absent in nonlinear CA. 

Further problematic features can be expected, when we consider the behaviour of nonlinear 
terms under the map relating the discrete description of CA and its continuum counterpart, 
employing Shannon’s Sampling Theorem as in Sec. 3., to which we turn shortly. 

4-2.1. A summary of properties of sinus cardinalis We will make use of several results 
concerning the sinus cardinalis function, sinc(x) := sin(x)/x, which is employed in the 
reconstruction formula, Eq. f|12D . and will be needed in the following. 

Let us define s n (t ) := sinc[ 7 r(t/t — n)], which leads to the Fourier transform: 

/ OO 

dt e~ lut s n (t) = lOfn/l — u)9(ir/l + u)e~' lU}ln , (19) 

-OO 

where 9 denotes the Heaviside step function. Thus, the function s n is bandlimited. Furthermore, 
it presents a “nascent” Dirac delta function, which has the properties: 

/ OO 

dt s n (t) = 1 , ( 20 ) 

-00 

/ OO 

dt s n (t)f(t) = 

-OO 


lim l 

l->0 


/( 0 ) • 


( 21 ) 


These results can be recovered with the help of the Fourier transform of s n , Eq. (1191) . assuming 
that / has a well-behaved Fourier transform. Employing the inverse Fourier transformation of 
Eq. (H9]), we obtain the orthogonality relation: 

/ OO 

d t s m (t)s n (t) = 5 mn . (22) 

-OO 

Finally, as an example, we perform a summation by applying the Sampling Theorem , Eq. (1121) : 

Y s n (ml - t')s m {t) = s n (t - t') , (23) 

mE Z 

noting that the factor of s n under the sum can be interpreted as the function on the right-hand 
side sampled at the times ml; all functions here have identical bandwidth, i.e. ui max = vr/ 1. 

4-2.2. Nonlinearity and nonlocality (in time) We now are ready to elaborate consequences 
of nonlinearity with respect to mapping the discrete equations of motion to the continuum 
desription thereof via application of Shannon’s Sampling Theorem. — Suppose that the discrete 
analogue of the Schrodinger equation, Eq. ©, includes a genuinely nonlinear term , e.g .: 

C = -iHapilfi + M a gry + V’nXC 7 + #0 > ( 24 ) 

keeping f n = 1 and where the coefficients are real and totally symmetric in the indices. 

A corresponding potential can be incorporated into the action and the additional nonlinearity 
in Eq. (l24l) follows by applying a suitably generalized variational derivative, as discussed in 
Sec. 4.1. We are not interested in the physical (ir)relevance of the present example of a cubic 
potential, but would like to see what happens with the nonlinear terms, e.g. M Q( g 7 V’n' i / ; n; when 
the Sampling Theorem is applied to Eq. (124[) . similarly as before in Sec. 3. 

Suppressing irrelevant greek indices, we introduce ip n =: t n ) and V’nV’n =: ^(2 )i)n)- 

Through the reconstruction formula (fT2l) the discrete time values V’(^n) and ^( 2 ){t n ) are replaced 
by continuous time functions ip(t) and ^( 2 )(f), respectively, and we would like to make the 
relation between the latter functions explicit. 

First of all, employing the orthogonality relation (1221) . we invert the reconstruction formula: 

/ OO 

d t s n (t)i/j(t) . (25) 

-OO 


Which gives us simply: 

/ OO p OO 

d f s n (t , )V’(^) / df" s n {t")il)(t") . (26) 

-oo J —OO 


Applying now the reconstruction formula to "0(2)(^n)> we obtain indeed a nonlinear relation 
between ^( 2 )(t) an d 


V , (2) (t)= l 2 



(27) 


where we interchanged summation and integrations. Similarly as in Eq. (1231) . we can do the sum: 
The function S( 2 ) n (t',t") := s n (t')s n (t") = sinc[7r(ni — t')/l]sinc[n(nl — t")/l\ is of the bandlimited 
kind and sampled here at the times nl; it is reconstructed by the summation including the 
factor s n (t), in accordance with the Sampling Theorem. However, the bandwidths need to be 


considered carefully. By Fourier transformation, one verifies that S( 2 )n has a doubled bandwidth, 

Umax = 27r/Z = 2u max , as compared to s n , which one would guess. This is implemented by 
writing all appearances of l in terms of 1/2 and by applying Eq. (1121) to yield the sum: 


S n(t)S( 2 )n(t' ,t") 

nGZ 


E r nl/2—t/2-\ 

smc F 1/2 \ 

nez 

• f Tll/2 — t* /2 "I • f 

• Sinc|7T ' t j 2 f Jsmc[7T 
sinc[7r(t — £')/(] sine [7r(t - 


nl/2—t"/2-\ 

1/2 J 


(28) 


Thus, we obtain: 

ip (2 = dt 7 sinc[7r(t — , (29) 

which expresses V’( 2 ) i n terms of ip. - In the limit of vanishing discreteness scale, we obtain a 
simple quadratic term : 

lirn -0 (2) (t) = {ip{t)) 2 , (30) 

with the help of Eq. (12TT) . This presents, of course, the expected result. It is local in time. 

However, we observe that the quadratic term on the right-hand side of Eq. (1291) consists in 
factors which are nonlocal in time: the function ip is integrated over all times, weighted by the 
oscillating and slowly decaying sine kernel. Inserting this result (and corresponding additional 
terms) into the continuous time version of the discrete analogue of a nonlinear Schrodinger 
equation, Eq. (1241) . changes the character of this equation profoundly: it is not anymore a 
consistent CA updating rule\ - This should be contrasted with the left-hand side of Eq. (fT3l) . 
which is nonlocal as well. However, this nonlocality is rather mild and refers to two neighbouring 
instants in such a way that the linear equation can be solved forwards (or backwards) in time 
step by step, recalling the discussion of initial conditions following Eq. (fT3l) in Sec. 3. With 
the nonlocality here, due to an anharmonic potential, inserted on the right-hand side of the 
continuous time equation, this fails. 

One may also consider that application of the Sampling Theorem to a nonlinear finite 
difference equation, such as Eq. (l24l) . cannot lead to a consistent continuous dynamics, since 
the resulting linear and nonlinear terms have different bandwidths (unless an additional cut-off 
on nonlinear terms is introduced by hand). Which holds for any kind of polynomial nonlinearity. 

Thus, we expect nonlocality to be a problem for any continuous description based on a form 
of sampling theory of an underlying discrete CA dynamics, unless the CA updating rules are 
linear in the dynamical variables (as in Sec. 2.). We anticipate this to be the case also if space 
is discrete, besides time, a situation which can be studied along similar lines [221 ]. 

This leaves us with a speculative question: Could it be that unitary linear evolution in 
continuous time — which appears to hold universally in QM (leaving aside measurement 
processes) — is dictated by a local perspective on more general, possibly nonlinear underlying 
CA dynamics? In short: Does locality filter for linearity? 


5. Conclusions 

We have briefly reviewed, in Sec. 2., the description of a class of deterministic discrete cellular 
automata based on an action principle PI El SI 1201. In particular, we have recalled how this can 
be mapped with the help of Shannon’s sampling theory HSUS] on a continuous time picture, 
which resembles in many respects the description of nonrelativistic quantum mechanical objects. 

In Sec. 3., we have pointed out the relation between the discrete CA updating rules, which are 
closely analogous to Hamilton’s equations of motion in mechanics, and a modified Schrodinger 





equation, which includes additional terms due to the finite discreteness scale l characterizing 
the CA (and leads to a modified dispersion relation, with energy bounded below and above). 
This extends to a one-to-one correspondence between the associated conservation laws, between 
continuous unitary symmetries and their discrete counterparts. 

Presently, in Sec. 4., we have paid particular attention to a generalization which incorporates 
polynomial nonlinearities into the action and equations of motion in a consistent way. 

This is motivated by earlier indications that only linear CA equations of motion could be 
consistent mm - which appeared a surprising result in view of the essential linearity of QM 
m- It is a fact that too naive implementation of genuine nonlinearities tends to produce 
overdetermined CA updating rules. 

We have found here that introducing a consistent nonlinear generalization of the discrete 
CA dynamics leads to an inconsistent nonlocality (in time) of the corresponding continuous 
time description obtained with the help of sampling theory. Furthermore, we have argued that 
nonlinearity has a detrimental effect (through Poisson brackets based on discrete variational 
derivatives) on the algebra of observables, as far as it compares with QM |21j . Both would 
severely spoil any attempt to construct models in the class of Hamiltonian CA which have QM 
features emerging at large scales (i.e. discreteness scale l —>• 0) from underlying CA dynamics. 

We conclude that nonlinearity must be excluded from the kind of CA model building studied 
here, in order to maintain locality in the continuous description to the extent possible, when 
ordinary derivatives are replaced by finite differences beneath. 

Yet one may wonder about effects of nonlinear CA processes, if they influence only some 
dynamical variables. Which could be different ones than those commonly described by a linear 
and local evolution in quantum theory. Is there room for stochastic phenomena that manifest 
themselves in QM measurements? 

In any case, in order to turn our observations of surprising connections between the properties 
of cellular automata and the quantum mechanical features of more familiar physical objects into a 
theory, as proposed by G. ’tHooft [1J, several immediate problems call for attention. We should 
understand how composite systems fare in this context. Which is a prerequisite to analyze 
the CA analogue of QM measurement processes and, more generally, the role of superposition 
principle and entanglement. Above all, we need to understand how aspects of relativity and of 
the physics of spacetime come into play here. 
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